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PsJ , Abstract. We study a dark matter halo model from two points of view: the "stellar polytrope" (SP) 

model coming from Tsallis' thermodynamics, and the one coming from the Navarro-Frenk- White 
■"nJ" I (NFW) paradigm. We make an appropriate comparison between both halo models and analyzing the 

^^ , relations between the global physical parameters of observed galactic disks, coming from a sample 

of actual galaxies, with the ones of the unobserved dark matter halos, we conclude that the SP model 
is favored over the NFW model in such a comparison. 
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INTRODUCTION 



An alternative formalism to the micro-canonical ensemble treatment, that allows non- 
extensive forms for entropy and energy under simplified assumptions, has been devel- 
oped by Tsallis [1] and applied to self-gravitating systems [2, 3, 4] under the assump- 
cd I tion of a kinetic theory treatment and a mean field approximation. As opposed to the 

Maxwell-Boltzmann distribution that follows as the equilibrium state associated with 
the usual Boltzmann-Gibbs entropy functional, the Tsallis' functional yields as equilib- 
rium state the "stellar polytrope" (SP), characterized by a polytropic equation of state 
with index «. 

On the other hand, high precision N-body numerical simulations based on Cold Dark 
Matter (CDM) models, perhaps the most powerful method available for understanding 
gravitational clustering, lead to the famous results of Navarro, Frenk and White (NFW 
model) [5] that predicts density and velocity profiles which are roughly consistent with 
observations, however, it also predicts a cuspy behavior at the center of galaxies that 
is not observed in most of the rotation curves of dwarf and LSB galaxies [6, 7, 8, 9, 
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10, 11]. The significance of this discrepancy with observations is still under dispute, 
nevertheless, the NFW model of collision-less WIMPs remains as a viable model 
to account for DM in galactic halos, provided there is a mechanism to explain the 
discrepancies of this model with observations in the center of galaxies. 

Since gravity is a long-range interaction and virialized self-gravitating systems are 
characterized by non-extensive forms of entropy and energy, it is reasonable to expect 
that the final configurations of halo structure predicted by N-body simulations must be, 
somehow, related with states of relaxation associated with non-extensive formulations 
of Statistical Mechanics; therefore, a comparison between the SP and NFW models is 
both, possible and interesting. The main purpose of our analysis is to make a dynamical 
analysis of two halo models, one based on the NFW paradigm, and other based on the 
SPs derived from Tsallis' non-extensive thermodynamics, compare them and test both 
with observational results coming from a sample of disk galaxies. 

SP AND NFW HALO MODELS 

For a face space given by (r, p), the kinetic theory entropy functional associated with 
Tsallis' formalism is [2, 3], and [4] 
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where / is the distribution function and c/ > 1 is a real number. In the limit q ^ \, 
the functional (1) leads to the usual Boltzmann-Gibbs functional, corresponding to 
the isothermal sphere. The condition 5Sq = Q leads to the distribution function that 
corresponds to the SP model characterized by the equation of state p = Knp^^^'", 
where Kn is a function of the polytropic index n, and can be expressed in terms of 

the central parameters: Kn = —fn- The polytropic index, n, is related to the Tsallis' 

Pc ' 

parameter q> \ by: n = | + -^. Inserting the equation of state into Poisson's equation, 
the Lane-Emden equation [12] is obtained, based on it we can obtain density, mass and 
velocity profiles for the SP model. 

NFW numerical simulations yield the following expression for the density profile of 
virialized galactic halo structures [5, 13]: 

^Po (2) 



Where: 5o = 3[in(i+coK"o/(i+co)] ' Po = Pcnt^o/^^ = 253.8/i2^o 0, y = coj;,, and Qq is 
the ratio of the total density to the critical density of the Universe. Using equation (2) its 
easy to obtain mass and velocity profiles. 



COMPARISON OF SP AND NFW MODELS 

In order to compare both halo models, it is important to make various physically mo- 
tivated assumptions. First, we want both models to describe a halo of the same scale, 
which means same virial mass M,,,. Secondly, both models must have the same maxi- 
mal value for the rotation velocity. This is a plausible assumption, as it is based on the 
TuUy-Fisher relation [14], a very well established result that has been tested success- 
fully for galactic systems, showing a strong correlation between the total luminosity of 
a galaxy and its maximal rotation velocity. Our third assumption is that the polytropic 
and NFW halos, complying with the previous requirements, also have the same total 
energy evaluated at the cut-off scale r = r,„. The main justification for this assumption 
follows from the fact that the total energy is a fixed quantity in the collapse and sub- 
sequent virialized equilibrium of dark matter halos [15]. Since the SP model has three 
free parameters (contrary to only one parameter of the NFW model), and we have se- 
lected three comparison criteria, we have a mathematically closed problem once M,„ is 
specified. 

Following the guidelines described above, we proceed to compare NFW and poly- 
tropic halos for M„, ranging from 10^*^ up to 10^^ solar masses. From this comparison 
we find the "best-fit" values for the free parameters of the SP model. The results are 
displayed explicitly in table 1 . The comparison between both models in velocity profiles 
is shown in the left panel of figure 1 . SP and NFW models have both the same virial 
mass, M,i, = lO^^M©. For other values of M,,, the velocity profiles are qualitatively simi- 
lar to the one displayed in figure 1 (left panel). The detailed description and results of the 
method of comparison presented above will appear in an article that is being prepared 
[16]. 

So far we have analyzed only the global structure of the dark matter halo without 
considering the effects of the luminous galaxy within. If one wishes to test a given model 
with observational results, it is necessary to add the galactic baryonic disk as a dynamical 
component of the model. In order to do so we followed the method described in [13, 
17]. Then to compare both models with observational results we used the prescription 
presented in [18]. 

Using such a prescription, we may define the ratio of maximum disk velocity to 
maximum total velocity, V(/,m/K.m, which is a global quantity that can be directly 
compared with theoretical predictions. This ratio is not defined at a given radius, but 
it can be related to the total mass to disk mass ratio Mt/Mj, defined at an specific radius. 
In particular at radius r„ where the total rotational curve has its maximum we have [18]: 



Md r^ \Vc 
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The use of Vd^m/Vcm instead of Mt/Mj is suitable because it can be obtained directly 
from observational parameters, without the assumptions needed to calculate Mt/Mj. 

One of the principal results obtained in the work [18] is that the ratio Vd.m/yc,m 
correlates principally with the disk surface density E^ of galaxies. Therefore we will 
use this result to compare the NFW and the SP models with the observational results 
coming from the sample. We will take M,,, = 1 x IO^^Mq, which is a characteristic value 



TABLE 1. Parameters characterizing the polytropes while being compared to NFW halos 

logio(M„,/M0) p,[Mq/^c^] C7,.[Km/s] n q Kn v^ax [Km/s] r.Jkpc] 
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4.82 


1.30 
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FIGURE 1. Left panel: velocity profiles for the NFW halo (solid line) and its compared SP one 
(dashed curve). Right panel: luminous to total dark matter content vs central surface density; open circles 
correspond to observational data, the dashed line represents the NFW model and the solid one to SP 



model. Both models have a virial mass of M,, = 10 M( 



for the mass of dark matter halos. For simplicity, we will assume that the baryonic mass 
fraction of the disk, /j, has the same value for all disks. 

In the right panel of figure 1 it is shown the ratio V^.m/K,™ vs LogLd for the sample 
of galaxies used (open circles). An almost linear trend can be seen between this quanti- 
ties. HSB (high surface brightness) galaxies (corresponding approximately to values of 
LogLd greater than 2.5) have greater values of Vd,ra/K-,m than LSB (low surface bright- 
ness) galaxies. This means that the luminous matter content is greater for HSB than for 
LSB galaxies. The shown picture is consistent with a well known result: LSB galaxies 
are dark matter dominated systems within the optical radius. 

The value of the graphic presented is that it allow us to bound statistically the possible 
values of the V^.m/K,™ ratio that galaxies with a given surface density can have. As 
was proved by [18], the size of this range of values (associated with dispersion on the 
graphic) is mainly due different virial mass that galaxies with the same E^ have. 

Right panel of figure 1 also show clearly that NFW models can not reproduce at 
satisfaction the results obtained for the compiled sample without introducing unrealistic 
values for the virial mass. This is one of the results that lead us to the possibility of 
seeking an alternative to the NFW paradigm. The curves shown in figure 1 (right panel) 
represent both models with average values for their respective parameter; it's clear from 
the figure that the SP model follows in better agreement the average behavior of the 
observational sample than the NFW model. 



CONCLUSIONS 

Motivated by the fact that SPs are the equilibrium state in Tsallis' non-extensive entropy 
formalism, we have found the structural parameters of those SPs that allows us to 
compare them with NFW halos of virial masses in the range 10^*^ < M,,JMq < 10^^; 
the results are displayed in Table 1. It is shown in the left panel of figure 1 that the 
velocity profile of the SP model is much less steep in the same region than that of 
the NFW halo. These features are consistent with the fact that NFW profiles predict 
more dark matter mass concentration than what is actually observed in a large sample of 
galaxies [19, 8, 18]. 

We have also shown that the SP model is favored over the NFW model regarding 
the dark matter content in disk galaxies (within the optical radius) which is shown by 
the average behavior of the observational sample in figure 1 (right panel). These results 
show that the NFW halo model can be enhanced with the use of alternative paradigms 
in Statistical Mechanics, which seems to solve a recurrent item which throws a shadow 
in such an excellent description as the NFW model is. 
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